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HIGHLIGHTS 


►  Modeling  of  thermal  and  capacity  fade  effects  for  a  LixC6-LiyMn204  cell. 

►  Capacity  fade  effects  at  anode  and  cathode  were  found  to  be  significant  during  cycling. 

►  SEI  formation  operates  under  steady-state  conditions  during  battery  cycling. 

►  Degradation  of  the  cell  is  more  prominent  during  the  dissolution  of  LiyMn204  particles. 
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Li-ion  batteries  are  the  most  promising  technology  for  use  in  electric  vehicles  in  the  near  future,  and  as 
such  it  is  critical  to  understand  their  performance  at  both  beginning  of  life  (BOL)  and  end  of  life  (EOL).  In 
this  work  different  thermal  and  capacity  fade  effects  (e.g.  SEI  formation,  dissolution  of  LiyMn204  parti¬ 
cles)  are  modeled  to  account  comprehensively  for  the  behavior  of  a  LixC6-LiyMn204  cell.  The  comparison 
between  baseline  and  complex  models  is  systematically  used  to  analyze  individual  contributions  and 
perform  a  deeper  evaluation  of  the  variables  affecting  the  capacity  fade  with  thermal  inputs  during 
typical  cycle  life  tests.  Some  modifications  in  the  original  model  are  proposed  to  better  describe  the 
behavior  of  the  cell  and  speed  up  the  calculations. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Interest  in  modeling  dual-lithium  insertion  cells  has  gained 
popularity  to  tackle  industrial  demands  because  of  its  efficiency  in 
testing  Li-ion  batteries.  Li-ion  batteries  are  the  most  promising 
technology  for  use  in  electric  vehicles  in  the  near  future.  Within  the 
automotive  environment  the  batteries  will  experience  demanding 
thermal  condition  (both  high  and  low  temperature),  large  state  of 
charge  swings,  and  high  power  draws.  All  of  these  harsh  operating 
conditions  have  the  potential  to  degrade  both  the  power  and 
capacity  performance  of  the  cells.  Therefore,  it  is  critical  to  under¬ 
stand  their  performance  at  beginning  of  life  (BOL)  and  end  of  life 
(EOL).  Modeling  is  also  a  powerful  method  to  analyze  failure 
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mechanisms  (e.g.  capacity  fade)  and  thermal  runaway  during  long 
cycle  life  tests  [1],  Arora  et  al.  [2]  reviewed  the  capacity  loss 
mechanisms  associated  with  phase  changes  in  the  inserted  elec¬ 
trode  materials,  electrolyte  decomposition,  active  material  disso¬ 
lution,  passive  film  formation,  etc.  Although  all  these  phenomena 
could  be  responsible  for  the  capacity  fade  in  batteries;  the 
consumption  of  lithium  at  the  anode  during  cell  charge  [2-4]  and 
cathode  dissolution  (e.g.  LiyMn204  spinels)  [4]  have  been 
acknowledged  as  the  most  dominant  effects.  In  this  direction, 
Spotnitz  used  modeling  to  propose  that  two  major  processes 
dominate  the  capacity  loss,  namely,  impedance  growth  at  both 
electrodes  and  loss  of  active  lithium  at  anode  [3],  Christensen  and 
Newman  also  suggested  that  the  use  of  cutoff  potentials  may  not 
prevent  undesired  parasitic  reactions  (e.g.  Li  salt  and  solvent 
consumptions  at  the  anode,  dissolution  of  LiyMn204  particles) 
occurring  during  overcharge  [4],  Ning  et  al.  [5]  modeled  the 
capacity  fade  of  Li-ion  batteries  considering  a  solvent  reduction 
(EC)  reaction  occurring  during  the  charging  process.  This  model 
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entailed  a  loss  of  active  lithium  ions  and  the  increase  of  the  anodic 
resistance  due  to  the  growth  of  a  film  (Solid  Electrolyte  Interface  or 
SEI).  More  recently,  Christensen  and  Newman  developed  a  model 
based  on  the  chemistry  of  the  SEI  (e.g.  film  properties)  and  the 
different  transport  phenomena  involved  during  its  formation  [6], 
Although  this  model  describes  the  film-growth  characteristics,  its 
application  is  limited,  because  it  does  not  account  for  other  coupled 
phenomena  occurring  across  the  cell.  This  limitation  is  likely  due  to 
the  complexity  of  the  film  formation  and  the  computational  efforts 
demanded  to  solve  the  entire  problem. 

Temperature  effects  on  the  performance  and  safety  of  the 
batteries  are  critical,  especially  in  automotive  applications  where 
battery  packaging  and  environmental  conditions  may  lead  to 
extreme  temperatures.  However,  unlike  capacity  fade,  thermal 
effects  are  subjected  to  more  inputs  (e.g.  heat  generation  due  to 
reactions,  entropy  changes,  Joule  heating)  during  cycling.  Different 
thermal  models  considering  cell  geometry  and  battery  stacks  have 
been  proposed  in  literature  for  the  Li-ion  system  [7,8].  However, 
none  of  them  has  considered  the  coupling  between  the  thermal 
effects  and  capacity  fade  effect  occurring  at  anode  and  cathode. 
Thus,  one  of  the  objectives  of  this  work  is  to  analyze,  through 
a  model  development,  the  interactions  arising  due  to  the  coupling 
of  thermal  and  degradation  phenomena  during  cycle  life  tests.  The 
model  presented  in  this  work  is  originally  based  on  the  model 
validated  with  experimental  data  conducted  by  Doyle  et  al.  [9],  the 
thermal  effects  reported  by  Gu  and  Wang  [7]  and  Cai  and  White  [8], 
the  capacity  fade  effects  due  to  an  electrochemical  reduction  at 
anode  by  Ning  et  al.  [5],  and  the  dissolution  of  LiyMn204  particles  in 
the  cathode  side  by  Park  et  al.  [10].  Comparisons  between  simple 
and  complex  models  are  carried  out  to  analyze  individual  degra¬ 
dation  effects  during  cycle  life  tests  considering  thermal  balances. 
Some  modifications  are  proposed  in  the  model  to  better  represent 
the  behavior  of  the  cell  and  accelerate  the  convergence  of  the 
model.  A  forthcoming  work  will  be  focused  on  determining  the 
phenomena  sensitivity,  the  accuracy  of  the  kinetics  parameters  of 
the  model  and  validating  the  model  with  the  experimental 
performance  of  the  cell. 


2.  Modeling 

The  dual-lithium  insertion  cell  modeled  in  this  work  is 
composed  of  a  negative  electrode  current  collector,  negative  elec¬ 
trode  (LixC6),  separator,  positive  electrode  (LiyM^O/t)  and  positive 
electrode  current  collector  as  shown  in  Fig.  1.  The  electrolyte  is 
a  solution  of  LiPF6  in  a  nonaqueous  liquid  mixture  of  ethylene 
carbonate  (EC)  and  dimethyl  carbonate  (DMC)  [9].  The  schematic  in 
Fig.  1  also  shows  the  dimensions  of  the  cross-section  of  the  cell 
modeled  in  this  study  with  a  1-D  domain.  The  intercalation  and 
deintercalation  reactions  taking  place  at  the  electrode/electrolyte 
interface  can  be  written  as: 


discharge 

Liy_xMn204  +  xLi+  +xe~  <  LiyMn2C>4  Cathode  (1) 

charge 

discharge 

LixC6  -c  Li0C6  +  xLi1  +  xe~  Anode  (2) 


During  discharge,  lithium  ions  insert  into  solid  particles  at  the 
cathode,  and  de-insert  from  particles  at  the  anode.  The  opposite 
occurs  during  the  charging  processes.  The  diffusive  process  in  the 
solid  phase  and  the  electrolyte  is  fundamental  in  the  mechanisms 
of  discharge  and  charge  of  the  cell,  as  well  as  the  heat  generation 
and  heat  transfer  across  the  cell.  The  governing  equations  for  the 
mathematical  model  presented  in  this  work,  as  well  as  the  initial 
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Fig.  1.  Schematic  of  the  dual  Li-ion  cell  showing  the  cross-section  (x-axis)  modeled  in 
this  work  (not  scaled).  The  cell  consists  of  a  negative  electrode  current  collector, 
negative  electrode  (LixC6),  separator,  positive  electrode  (Lij,Mn204)  and  positive  elec¬ 
trode  current  collector.  The  electrolyte  is  a  solution  of  LiPF6  in  a  nonaqueous  liquid 
mixture  of  ethylene  carbonate  and  dimethyl  carbonate. 


and  boundary  conditions  required  to  solve  the  equations  are  shown 
in  Table  1.  The  mathematical  model  describes  the  diffusion  (e.g. 
porous  electrode  theory)  and  conduction  of  Li  ions  with  conser¬ 
vation  of  charge  (e.g.  Ohm’s  law)  in  the  solid  and  electrolyte  phases 
across  the  cell,  e.g.  anode,  separator  and  cathode.  Further  details  of 
the  material  and  charge  balances  have  been  described  elsewhere 
[9],  The  model  also  incorporates  energy  balances  to  describe  the 
thermal  distributions  across  the  cell  [7,8],  the  equations  to  describe 
capacity  fade  effects  due  to  the  loss  of  active  lithium  ions  by  elec¬ 
trochemical  reduction  [5]  and  volume  changes  due  to  the  dissolu¬ 
tion  of  LiyMn204  particles  at  the  cathode  side  [10],  Note  that 
a  modification  was  applied  to  the  model  to  transform  the  discrete 
expressions  reported  in  Ref.  [5]  for  the  lithium  concentration  in  the 
anode  and  resistance  of  SEI  formation  into  derivative  terms.  These 
changes  remarkably  improved  the  convergence  of  the  model. 
Table  2  details  the  parameters  and  formulas  used  in  the  lithium-ion 
battery  model  reported  in  Table  1. 

Numerical  simulations  were  run  in  order  to  solve  the  model 
reported  in  Table  1.  Matlab™  2010b  and  COMSOL™  MULTIPHYSICS 
4.1  were  used  together  for  this  purpose.  Second-order  Lagrange 
quadratic  polynomials  were  used  as  shape  functions,  along  with  an 
interactive-smart  mesh  used  to  discretize  the  domain  using  the 
Finite  Element  Method.  An  Intel  Core  i7  central  processing  unit 
running  at  3.47  GHz  with  a  random  access  memory  of  8.00  GB  was 
used  to  carry  out  the  calculations. 

3.  Results 

The  first  set  of  simulations  performed  to  test  our  model  were 
comparisons  of  potential  vs.  capacity  of  the  cell  at  different 
discharge  rates  (plots  not  shown).  These  computations  were  run 
under  the  conditions  reported  by  Doyle  et  al.  (refer  to  Fig.  4  in  Ref. 
[9]).  Our  simulations  showed  that  the  model  was  capable  of 
describing  the  experimental  plots  shown  by  Doyle  et  al.,  and  only 
slight  deviations  were  observed  likely  related  with  the  choice  of  the 
numerical  method  used  to  solve  the  model.  For  example,  Doyle 
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et  al.  used  Finite  Differences  instead  of  Finite  Element  Method.  The 
intention  of  this  work  is  to  analyze  the  degradation  and  the  inter¬ 
actions  arising  due  to  the  coupling  of  capacity  fade  and  thermal 
effects  over  the  life  of  a  cell.  Therefore,  the  next  step  was  to  cycle 
the  battery  using  the  typical  cycling  profile  shown  in  Fig.  2.  A 
baseline  model  (i.e.  without  coupled  effects)  was  compared  with 
more  comprehensive  ones  to  analyze  systematically  in  the  effects 
of  the  following  three  combined  factors  on  cell  performance:  1) 
thermal  balance,  2)  thermal  balance  +  capacity  fade  at  the  anode, 
and  3)  thermal  balance  +  capacity  fade  at  the  cathode. 

3.1  Current  density  profile 

Fig.  2  shows  the  current  density  profile  (i.e.  single  cycle)  for  the 
charge-depleting  cycle  life  test  reported  by  the  Idaho  National 
Laboratory  in  the  battery  test  manual  for  plug-in  hybrid  electric 
vehicles  [11],  This  waveform  is  intended  to  allow  the  vehicle  to 
operate  in  hybrid  and  electric  modes  with  a  net  decrease  in  the 
battery  state  of  charge  (SOC)  and  with  an  observable  degradation  of 
the  battery  from  cycle  to  cycle.  Thus,  this  typical  cycle  fulfills  the 
requirements  to  test  the  model  that  has  been  developed  in  this 
work.  This  profile  was  used  to  run  all  the  simulations  reported  in 
this  work.  Note  that  for  batteries,  this  profile  assumes  that  the 
degradation  rate  is  a  function  of  the  high  current  densities  passing 
through  each  cell.  Then,  the  current  passing  through  each  cell  will 
be  proportional  to  the  total  cell  cross-sectional  area  (which  will 
depend  on  how  the  cell  is  engineered),  and  the  number  of  cells  in 
parallel.  Therefore,  for  this  study  it  is  recognized  that  the  degra¬ 
dation  rate  is  high  with  the  intention  of  enhancing  the  degradation 
effects  over  cycling  (e.g.  capacity  higher  than  50  C).  In  future  work, 
the  sensitivity  of  this  rate  and  the  current  passing  through  each  cell 
will  be  explored  in  more  details. 

3.2.  Thermal  contributions  without  capacity  fade  effects 

The  first  cycling  tests  using  the  profile  in  Fig.  2  were  computed 
to  analyze  the  thermal  effects  on  the  behavior  of  the  battery  cell. 
For  this  phase  of  the  work  different  heat  transfer  coefficients  (h)  as 
well  as  isothermal  and  adiabatic  conditions  were  used  in  the 
model.  These  results  are  shown  in  Fig.  3  for  50  cycles  (equivalent  to 
18,050  s  on  cell  operation).  No  capacity  fade  effects  were  consid¬ 
ered  at  this  stage  of  the  modeling  and  the  linear  dependence  of  the 
OCP  is  neglected,  dUn/dT  —  0  [7],  As  shown  in  Fig.  3,  the  highest 
temperatures  are  reached  when  the  cell  is  placed  under  adiabatic 
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cycle 

Fig.  3.  Computed  temperature  profiles  at  the  interface  formed  between  the  cathode 
and  the  current  collector  at  different  heat  transfer  coefficients  labeled  on  the  figure. 


conditions,  followed  by  when  the  process  is  performed  under  the 
lowest  heat  transfer  coefficient,  h  —  0.01  W  m-2  K-1.  These  effects 
are  not  surprising  since  the  cell  lies  under  isolated  conditions  and 
low  heat  exchange  with  the  environment,  respectively  in  each  one 
of  these  simulation  runs.  On  the  other  hand,  when  the  cell 
exchanges  heat  easily  with  the  environment  (fi>lW  m2  K-1),  it 
attains  a  steady  temperature  profile  over  cycling  since  there  is  not 
sufficient  heat  generated  to  cause  a  rise  in  the  temperature. 
Typically  h  is  assumed  to  be  1  W  m  2  I<  1  but  there  is  expected  to 
be  a  wide  variation  in  this  value  depending  on  how  the  cell  is 
manufactured.  Two  regions  are  of  interest  under  these  conditions. 
In  the  first  region,  during  the  first  236  s  of  each  cycle  (low  current 
zone  in  Fig.  2),  the  heat  produced  in  the  cell  is  mostly  exchanged 
with  the  surroundings.  The  opposite  behavior  occurs  in  the  second 
region  ( ~  last  236-268  s  of  each  cycle)  when  larger  discharge/ 
charge  current  densities  are  imposed  on  the  cell,  producing 
enough  heat  to  increase  the  temperature  within  the  cell.  This 
observation  is  repeated  through  each  cycle  and  a  steady  temper¬ 
ature  profile  is  maintained.  The  rises  and  drops  in  temperature 
shown  in  the  inset  of  Fig.  3  correspond  to  abrupt  and  intermediate- 
low  Li+  intercalation/deintercalation  processes,  respectively  (plots 
not  shown).  It  is  worth  highlighting  that  steeper  slopes  are 
attained  in  the  temperature  profiles  when  higher  currents  are 
applied  to  the  cell.  Apparently,  these  effects  would  multiply  when 
battery  stacks  are  cycled  over  the  long  term,  whence  thermal 
overheating  would  be  also  produced.  It  is  worth  mentioning  that 
these  simulations  do  yet  to  account  for  effects  that  may  be 
observed  near  the  current  collectors  where  higher  localized 
current  densities  may  be  observed  (e.g.  Joule  heating  effect).  The 
variation  of  the  temperature  across  the  cell  was  found  to  be 
insignificant  regardless  of  the  condition  utilized  in  the  cell.  This 
behavior  could  be  associated  with  the  fact  that  the  present  model 
considers  only  the  cross-section  of  the  cell.  Larger  effects  might  be 
expected  along  the  height  of  the  cell  where  Joule  heating  effects 
may  become  important  [7], 

3.3.  Capacity  fade  effects  due  to  anode  SEI  formation 

In  Fig.  4  anode  degradation  is  assumed  from  SEI  formation,  and 
it  is  assumed  that  the  SEI  continues  to  increase  with  each  cycle.  As 
such.  Fig.  4  shows  the  computed  potential  profiles  at  different 
exchange  current  densities  of  the  electrolyte  reduction  reaction 
occurring  during  the  charge  process  at  the  anode  (i.e.  different 
rates  for  SEI  formation).  This  reaction  is  used  to  account  for  the 
capacity  fade  effects  in  the  battery  due  to  the  growth  of  the  SEI  [5], 


a  continuous  formation  of  SEI  film  throughout  the  cell  life.  Different  exchange  current 
densities  accounting  for  SEI  formation  are  evaluated  using  a  heat  transfer  coefficient 
set  to  1  W  m  2  JT1.  Cycles  0-50. 


Thermal  effects  are  also  included  in  this  model.  As  observed  from 
these  plots,  the  Eceii  oscillates  at  the  two  highest  j'sei  values  (0.08 
and  0.0008  A  m-2),  probably  due  to  the  significant  increase  of  the 
reduction  rate  for  this  reaction  as  the  number  of  cycles  are 
increasing  (jsEi )•  This  process  causes  a  loss  of  active  lithium  and  an 
increase  of  impedance  ( RSei )  at  the  anode  (plots  not  shown),  as 
well  as  the  oscillation  of  the  charge  in  the  solid  phase  (refer  to  $i 
in  Table  1 )  and  a  minor  oscillation  of  the  electric  potential  in  the 
liquid  phase  ($2)-  On  the  other  hand,  when  the  capacity  fade  is  not 
considered  in  the  model  (i.e.  baseline  model)  or  low  contributions 
of  the  reduction  reaction  are  considered  (j'sei  =  0.8  x  \0  1  A  m~2), 
the  Eceii  profiles  resemble  those  obtained  experimentally  (plots 
not  shown).  and  $2  do  not  oscillate  for  these  former  cases. 
However,  the  simulations  carried  out  with  the  baseline  model, 
using  jsei  =  0.8  x  10-7  or  0.8  x  10-5  A  m-2  (not  shown)  produce 
virtually  the  same  Ecen  profiles.  This  behavior  is  associated  with 
the  low  j'sei  values  used  in  the  simulations,  which  are  not  sufficient 
to  generate  an  important  input  of  the  SEI  formation.  As  previously 
described,  larger  j'sei  values  lead  to  the  oscillation  of  Eceu,  which 
would  not  represent  a  realistic  description  of  the  behavior  of  the 
cell.  This  phenomenon  suggests  that  the  formation  of  the  SEI  at 
the  anode  would  entail  a  different  action.  Perhaps,  its  formation 
might  start  occurring  during  the  early  cycles  of  the  test,  and  reach 
a  steady-state  for  the  consecutive  cycles.  In  fact,  during  the 
manufacturing  of  Li-ion  batteries,  they  are  initially  cycled  for 
some  cycles  to  form  a  pseudo-steady-state  SEI,  and  thus  avoiding 
the  loss  of  Li+  ions  from  the  electrolyte  during  subsequent  oper¬ 
ation.  Under  the  scheme  in  Fig.  4,  a  steady  loss  of  active  lithium 
and  impedance  growth  would  be  also  implied  in  the  model.  In 
fact,  important  phenomena  associated  with  the  SEI  formation  are 
apparent  from  literature  [4],  Numerous  reactions  could  be 
involved  in  the  SEI  formation  when  the  anode  contacts  an  oxidant 
(e.g.  solvent,  air)  [6],  whence  the  film  composition  has  been 
difficult  to  determine  [2],  Nevertheless,  it  is  known  that  the  film  is 
to  some  extent  electronically  insulating  whereby  continuous 
solvent  and  salt  reduction  at  the  S£//solution  interface  is  hindered 
[6],  Thus,  stable  films  form  on  graphite  electrodes  during  the  first 
two  cycles  [12]  and  any  further  change  of  their  structure  occurs 
only  under  abusive  conditions  or  long  storage  [2],  Solvent  mole¬ 
cules  are  impeded  from  arriving  at  the  anode,  whereas  Li+  ions  are 
transferred  across  the  film  as  vacancies  or  interstitials  [6]. 
Accordingly,  the  formation  of  a  steady  film  after  the  first  two 
cycles  is  used  below  to  describe  the  capacity  fade  effects  during 
cycle  life  tests. 
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Fig.  5  shows  Iap  (profile  in  Fig.  2)  and  Eceu  profiles  for  cycles 
49—50  assuming  that  the  SEI  thickness  and  consequently  the 
impedance  growth  remain  unchanged  after  the  first  two  cycles 
according  to  information  provided  in  the  literature  [2,4,13].  It  is 
acknowledged  that  that  for  some  particular  anodes,  the  SEI 
continues  growing  but  its  effects  of  the  addition  growth  after  the 
initial  SEI  formation  can  be  considered  negligible  [14],  unless  it  g 
overcomes  a  major  damage  which  require  a  big  loss  of  Li+  ions  to  111 
repair  it.  Instead  of  the  discharge/charge  capacity  plots,  Eceu  profiles 
are  presented  since  cycles  containing  different  waveforms  were 
applied  in  this  study.  The  curves  in  Fig.  5  were  overstated  by  the  use 
of  large  values  of  j'se/  to  appreciate  the  effects  of  capacity  fade. 
However,  none  of  them  present  an  oscillatory  behavior  as  those 
shown  in  Fig.  4.  The  difference  is  evident  between  the  model  that 
does  not  consider  a  capacity  fade  input  and  the  one  that  accounts 
for  this  effect.  Note  that  the  baseline  model  produces  virtually 
similar  results  with  the  simulations  carried  out  at 
Jsei  -  0.8  x  10  7  A  m  2,  which  is  likely  associated  with  a  low  input 
of  this  exchange  current  density.  As  observed  in  the  discharge 
regions  of  the  profile,  the  Eceu  is  in  general  lower  for 
jsEi  —  0.8  x  1CT  A  m-2,  following  this  trend  the  value  of 
0.8  x  10-3  A  m  2.  This  suggests  that  higher  jsei  values  accelerate  the 
SEI  precipitation,  whereby  the  film  resistance  presented  to  lithiate/ 
delithiate  Li+  ions  increases  during  the  first  two  cycles.  Under  this 
context,  larger  j'se;  values  would  cause  a  bigger  capacity  fade  of  the 
cell.  Likewise,  Fig.  5  clearly  shows  the  prediction  of  the  accelerated 
degradation  imposed  to  the  cell  over  cycling.  Thus,  since  oscillatory 
behaviors  are  obtained  for  the  Eceu  when  the  SEI  resistance  is 
allowed  to  grow  (refer  to  Fig.  4),  the  impedance  of  the  film  should 
remain  steady  during  most  of  the  cycling. 

3.4.  Capacity  fade  effects  due  to  changes  in  the  volume  fraction  of 
the  cathode 

Considering  cathode  degradation,  Fig.  6  shows  two  different 
energy  activation  (Ea)  values  describing  the  volume  change  due  to 
dissolution  of  the  LiyMn2C>4  particles  at  the  cathode.  The  value  of 
72,840  J  mol-1  reported  for  this  process  was  selected  from  Ref.  [10], 
and  it  was  only  decreased  by  an  order  of  magnitude  to  sensitize  its 
contribution.  Fig.  6a  shows  a  clear  difference  in  Eceu  values  when 
two  different  Ea  values  are  used  in  the  simulations  during  the  first 
five  cycles.  Not  surprisingly,  lower  Eceu  values  are  obtained  when 
the  energy  activation  to  dissolve  LiyMn204  particles  is  lower,  indi¬ 
cating  that  the  cell  capacity  degrades  faster  when  the  energetic 
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Fig.  6.  Computed  cell  potential  for  analyzing  the  capacity  fade  at  cathode  assuming 
a  change  in  volume  fraction  of  this  electrode  due  to  the  dissolution  of  the  LiyMn204 
particles.  Two  different  energy  activation  values  are  evaluated  using  a  heat  transfer 
coefficient  set  to  h  =  1  W  m~2  K-1.  Capacity  fade  at  anode  is  disregarded.  Cycles  1-5 
and  40-45  are  shown  in  (a)  and  (b),  respectively. 

barrier  to  perform  this  reaction  is  decreased.  A  comparison  of  the 
Eceu  profiles  given  for  cycles  1  —5  (Fig.  6a)  and  40—45  (Fig.  6b)  shows 
that  a  larger  difference  between  the  two  Ea  values  is  observed  as  the 
cycling  of  the  cell  becomes  longer,  whereby  the  degradation  of  the 
cell  progressively  increases  with  cycling.  Variations  in  the  standard 
rate  constant  of  the  dissolution  of  LiyMn2C>4  particles  (ko)  did  not 
produce  differences  in  the  computed  Eceii.  This  might  be  due  to  the 
fact  that  k0  and  Ea  are  statistically  correlated  in  the  kdiss  expression 
(refer  to  Table  2),  and  k0  is  not  sensitive  to  the  model  under  the 


-50  g 
< 

-100 


cycle 

Fig.  5.  Computed  cell  potential  for  analyzing  the  capacity  fade  at  anode  assuming  the 
formation  of  SEI  film  only  occurs  during  the  first  two  cycles  and  then  remains  constant. 
Different  exchange  current  densities  accounting  for  SEf  formation  are  evaluated  using 
a  heat  transfer  coefficient  set  to  1  W  m2  K  Simulations  of  the  baseline  model  and 
Jsei  =  0.8  x  10  7  A  nrr2  produce  the  same  result  Cycles  49-50. 
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analyzed  conditions.  The  results  obtained  for  this  specific  degra¬ 
dation  mode  need  to  be  experimentally  validated  in  order  to 
determine  the  contributions  of  the  dissolution  of  the  LiyMn204 
particles  operating  in  real  Li-ion  batteries  (e.g.  any  potential 
correlation  involving  heat  generation).  Presumably,  other  impor¬ 
tant  factor  that  should  be  analyzed  is  the  migration  and  subse¬ 
quently  deposition  of  Mn2+  ions  in  the  anode,  which  could  also  lead 
to  the  capacity  fade  of  the  battery. 

Fig.  7  shows  the  volume  fraction  (V//Vfni)  changes  occurring  in 
the  cathode  due  to  the  dissolution  of  the  LiyMn204  particles 
when  an  energy  activation  value  of  72,840  J  mol-1  is  considered 
in  the  simulations.  As  observed  in  this  plot,  the  volume  fraction 
slightly  augments  linearly  as  cycling  progresses.  This  behavior 
shows  a  similar  trend  described  in  the  literature  to  a  battery  that 
is  stored  at  20  °C  for  seven  days,  and  larger  volume  fraction 
changes  are  observed  when  the  battery  is  stored  at  50  °C  during 
the  same  time  [10], 

4.  Conclusions 

A  comprehensive  model  has  been  presented  to  analyze  the 
interactions  between  thermal  and  capacity  fade  effects  during 
cycling  of  Li-ion  cells.  Heat  generation  due  to  reactions,  entropy 
changes  and  Joule  heating  are  considered  for  the  thermal  effects, 
whereas  the  capacity  fade  effects  are  accounted  for  by  the  loss  of 
active  lithium  ions  through  electrochemical  reduction  at  the  anode 
(SEI  formation),  and  the  volume  change  due  to  dissolution  of 
LiyMn204  particles  at  the  cathode.  Some  modifications  in  the 
capacity  fade  terms  were  carried  out  in  order  to  improve  the 
convergence  of  the  model  and  account  for  the  behavior  of  the  cell. 
The  baseline  model  presented  in  this  work  (e.g.  mass  and  charge 
balances)  agreed  with  the  behavior  presented  in  other  studies. 

Thermal  effects  analysis  revealed  that  the  cell  exchanged 
considerably  heat  for  h  >  1  W  m  2  K-1,  attaining  a  stable  temper¬ 
ature  profile  during  cycling.  On  the  other  hand,  adiabatic  conditions 
imposed  to  the  cell  produced  a  progressive  rise  in  temperature 
since  isolated  conditions  prevailed  in  the  system.  This  effect  was 
promoted  by  the  application  of  high  currents  in  the  cell.  The 
thermal  analysis  showed  the  need  to  incorporate  these  effects  in 
the  modeling  of  Li-ion  cells,  in  addition  to  the  integration  of  the 
heat  produced  by  different  cells  in  a  battery. 


Computed  potential  (Ece»)  profiles  were  used  to  discriminate  the 
capacity  fade  effects  due  to  SEI  formation.  Oscillatory  potential 
profiles  were  obtained  when  this  reaction  was  considered  tran¬ 
sient,  and  the  impedance  of  the  film  grew  over  cycling.  Otherwise, 
the  Ecen  profiles  resemble  the  behavior  obtained  experimentally 
when  these  inputs  were  considered  transient  during  the  first  two 
cycles,  and  steady  for  the  rest  of  the  test.  Results  from  the  evalua¬ 
tion  of  the  capacity  fade  at  cathode  showed  that  the  degradation  of 
the  cell  is  more  prominent  by  using  a  lower  energy  activation  value 
used  to  account  for  the  dissolution  of  LiyMn204  particles.  Since 
capacity  fade  effects  at  anode  and  cathode  were  found  to  be 
significant  in  the  degradation  of  the  cell  under  specific  conditions. 
It  is  needed  to  compare  the  model  predictions  with  experimental 
data  modulated  over  different  degradation  rates  to  be  able  to 
elucidate  the  importance  of  these  effects  under  those 
circumstances. 
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